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Abstract 

The cross section of W-boson pair-production is measured with the L3 detector 
at LEP. In a data sample corresponding to a total luminosity of 629.4 pb~^, collected 
at centre-of-mass energies ranging from 189 to 209 GeV, 9834 four-fermion events 
with W bosons decaying into hadrons or leptons are selected. The total cross section 
is measured with a precision of 1.4% and agrees with the Standard Model expec- 
tation. Assuming charged-lepton universality, the branching fraction for hadronic 
W-boson decays is measured to be: Br(W hadrons) = 67.50 ± 0.42 (stat.) ± 
0.30 (syst.) %, in agreement with the Standard Model. Differential cross sections 
as a function of the W~ production angle are also measured for the semi-leptonic 
channels qqeu and qq/iz/. 
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1 Introduction 



From 1996 until the year 2000, the centre-of-mass energy, y^, of the LEP e"''e~ colhder at 
CERN was increased in several steps from 161 GeV to 209 GeV. These energies, being above 
the kinematic threshold of W-boson pair production, allow detailed studies of this process. 

To lowest order within the Standard Model [1], three chargcd-current Feynman diagrams, 
shown in Figure 1 and referred to as CC03 [2-4], yield four-fermion final states via W-boson 
pair production: t-channel v exchange and s-channcl 7 and Z-boson exchange. W bosons decay 
into a quark- antiquark pair or a lepton-antilepton pair, denoted here as qq, {i = e, /x, r) or, 
in general, ff for both W+ and W~ decays. This Letter describes measurements of all four- 
fermion final states ipiu, qqiu and qqqq mediated by W-boson pair production. The presence 
of additional photons in the final state is not excluded. Contributions to the production of four- 
fermion final states arising from other neutral- or charged-current Feynman diagrams are small. 
At the current level of statistical accuracy, interference terms are sizable for the 151 charged- 
and neutral-current diagrams contributing to the ii/iu final states, for the 20 charged-current 
diagrams contributing to the qqeu final state and for the 214 charged- and neutral-current 
diagrams contributing to the qqqq final state [2-4] . 

It is conventional to quote results for the CC03 subset of diagrams, including the effect of 
initial-state radiation. As four-fermion states produced by all diagrams are measured, a suitable 
reweighting technique, described in the following, is applied to extract these results. 

W-boson decay branching fractions and the total W-boson pair-production cross section 
are determined with improved precision as compared to earlier L3 measurements at ^/s — 
161 — 189 GeV [5-8]. Comparable results were reported by other LEP experiments [9]. 

The differential cross sections for the qqeu and qqnv final states, as a function of the W~ 
production angle with respect to the direction of the incoming electron, are also derived. 

2 Data and Monte Carlo Samples 

The results presented in this Letter are based on the full luminosity collected by the L3 de- 
tector [10] during the high-energy runs of the LEP collider. The measurement of the total 
luminosity, C, follows the procedure described in Reference 11. 

The data collected at ^/s — 192 — 209 GeV are analysed in seven ^/s bins, as detailed 
in Table 1. The corresponding centre-of-mass energies are known with a precision of about 
50 MeV [12]. Results based on data collected at ^/s = 189 GeV were already published [8] but 
are reanalysed here since improved Monte Carlo programs arc now available for signal simu- 
lation: KandY [13] and RacoonWW [14]. The KandY generator combines the four-fermion 
generator KORALW [15] with the 0{a) corrections to W-boson pair production as imple- 
mented in the YFSWW3 [16] program. These corrections are calculated using the leading-pole 
approximation [17]. The RacoonWW Monte Carlo program implements such corrections in the 
double-pole approximation with similar accuracy and it is used for the estimation of systematic 
uncertainties. 

All KandY Monte Carlo samples used in this analysis are generated using the full set 
of Feynman diagrams contributing to a specific four-fermion final state. The KandY program 
provides the matrix elements on an event-by-event basis for different contributions including, for 
example, the CC03 subset of diagrams or the 0{a) corrections. This feature is extensively used 
in the following, both to derive quantities at CC03 level and for the assessment of systematic 
uncertainties. For example, the CC03-level efficiencies are calculated by reweighting every 
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event with the factor wccos = |7W(CC03)|7|A^(4f)p, where 7W(CC03) and A^(4f) are the 
matrix elements of the CC03 subset of diagrams and of the full set, respectively. The same 
events, reweighted by the factor 1 — WccoS) describe the remaining four-fermion contribution 
not arising from W-boson pair production. In the following they are considered as background. 
As a cross-check, selection efficiencies are also derived using the EXCALIBUR [18] four-fermion 
generator. 

The following Monte Carlo event generators are used to simulate the background processes: 
KK2f [19], PYTHIA [20], BHAGENE3 [21] and BHWIDE [22] for fermion-pair production, 
denoted as e+e" ^ ff(7); TEEGG [23] for radiative e+e" e+e--f{-f) events; DIAG36 [24] 
and LEP4F [25] for two-photon collisions with lepton-pair final states and PHOJET [26] for 
two-photon collisions with hadronic final states. 

Quark fragmentation and hadronisation processes are simulated using PYTHIA. Its param- 
eters are tuned to describe hadronic Z decays at y/s = 91 GeV [27]. A dedicated parameter set, 
derived from a light-quark Z-decay data sample, is used for the W-boson pair-production sim- 
ulations. Bose- Einstein correlations between hadrons from W decays are simulated using the 
LUBOEI BE32 model [28], with Bose- Einstein correlations only between hadrons originating 
from the same W boson, as supported by our study [29]. 

The response of the L3 detector is modelled with the GEANT [30] detector simulation pro- 
gram which includes the effects of energy loss, multiple scattering and showering in the detector 
material. Hadronic showers are simulated with the GHEISHA [31] program. Time-dependent 
detector inefficiencies, as monitored during data taking, are included in the simulations. 

3 Four-Fermion Event Selection 

The selections of four-fermion final states are designed to mimimise the uncertainty on the cross 
section of each channel. They are chosen to be mutually exclusive, by using complementary 
cuts, in order to avoid double counting of events. 

Electrons are identified as energy depositions in the BGO electromagnetic calorimeter having 
an electromagnetic shower shape and matching in azimuth a track reconstructed in the central 
tracking chamber. Muons are identified as tracks reconstructed in the muon chambers, which 
point back to the interaction vertex. Tracks which match a minimum-ionising-particle signature 
in the calorimeters are also retained as muon candidates and denoted as MIPs. Jets arising from 
hadronic tau decays are reconstructed using a jet-clustering algorithm in a cone of 15° half- 
opening angle [32] . The momentum of the neutrino in qqiv events is identified with the missing 
momentum vector of the event. Hadronic jets corresponding to quarks are reconstructed using 
the Durham jet algorithm [33]. In the e~'"e~ qqiu selections, the hadronic jets are formed 
from energy depositions and tracks not belonging to the reconstructed lepton. 

Efficiencies are evaluated for each ^/s point in the form of 10 by 10 matrices relating events at 
CC03 level to those at reconstruction level. An example is given in Table 2 for ^/s — 206.5 GeV. 
Selection efficiencies at other centre-of-mass energies are only marginally different. 

The number of selected events and background contributions are detailed in Table 3. A 
more detailed description of all selections is given below. 
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3.1 The iuiu Selection 



The event selection for the process e+e~ — > Ivlv requires two charged leptons and missing 
energy due to the neutrinos. The selection depends on whether the event contains zero, one or 
two identified electrons or muons, referred to as jet-jet, lepton-jet and lepton-lepton classes. For 
the lepton-jet and jet-jet classes, only the most energetic jets are retained as tau candidates. 
Electrons, muons and jets from hadronic tau decays are identified within the polar angular 
range | cos6'| < 0.96, where 6 is the lepton angle with respect to the beam direction. For events 
with one or two electrons, one electron is required to satisfy | cos^| < 0.92 in order to reduce 
the background from Bhabha scattering. For the jet-jet class, the two most energetic jets must 
also satisfy | cos^| < 0.92. 

The acoplanarity, defined as the complement of the angle between the directions of the two 
leptons in the plane transverse to the beam direction, must be greater than 8 degrees for the 
lepton-lepton and lepton-jet classes and 14 degrees for the jet-jet class. These criteria suppress 
the dominating backgrounds from lepton-pair production and cosmic rays. The leptons must 
have a signal in the scintillator time-of-flight counters compatible with the beam crossing. The 
total momentum transverse to the beam direction, P^, must be greater than 8 GeV. 

Events belonging to the lepton-lepton class are selected by requiring an energy of at least 
25 GeV for the more energetic lepton and 5 GeV for the less energetic one. For the lepton-jet 
class, the energies of the lepton and of the jet must exceed 20 GcV and 8 GeV, respectively. 
For the jet-jet class, the energies of the most energetic and second most energetic jets must be 
greater than 20 GeV and 6 GeV, respectively. 

The selected sample has a purity of 72% at -y/s = 206.5 GeV. The remaining background 
is dominated by lepton production in two- photon collisions (50%) and lepton-pair production 
(24%). The distributions of the acoplanarity and of the missing momentum transverse to the 
beam direction for the lepton-lepton class are shown in Figure 2. 

3.2 The qqei/ Selection 

The event selection for the process e+e" — > qqei/ requires an identified electron of at least 
20 GeV, high particle-multiplicity and large missing momentum. 

The missing momentum direction must point well inside the detector with a polar angle, 
6'miss, such that | cos6'miss| < 0.95. The reconstructed jet-jet and lepton-neutrino masses, referred 
to as Mjj and Mej,, must be greater than 45 GeV and 63 GeV, respectively. The latter cut is 
used to discriminate between e"'"e~ — > qqez/ and e"'"e~ — > qqrz/ events with r — >■ evv. To 
further suppress the dominant background from the e^e~ qq(7) process, which is planar, 
the directions of the electron and of the two jets are required to subtend a solid angle of less 
than 5.3 sr. 

The purity of the selection is 98% at ^/s — 206.5 GeV. The accepted background not 
originating from W-boson pair production is dominated by e"'"e~ — > qqez/ final states (71%) and 
e"'"e~ qq(7) events (29%). The distributions of the energy of the electron and of | cos(^niiss)| 
are shown in Figure 2. 

3.3 The qqpLU Selection 

The event selection for the process e"'"e~ — > qq/xz/ requires high particle-multiplicity, an identified 
muon or a MIP, and large missing momentum. 
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The jet-jet mass must satisfy 25 GeV < Mjj < 125 GeV for events with identified muons and 
50 GeV < Mjj < 98 GeV for events with MIPs. The muon-neutrino reconstructed mass, M^,^, 
is used as a discriminant against e+e^ qqru events with r ^iVt^^h- Its value is required to 
exceed 53 GeV. This cut is not apphed for events containing MIPs. 

The discrimination against e'^e" — > qqrv events is further enhanced by requiring the variable 
-P* — \Pix\ — 10 GeV(cos^* + 1), where is the momentum of the muon and d* is the decay 
angle of the muon in the reconstructed W-boson rest frame, to satisfy P* > 18.5 GeV. This 
requirement is loosened to P* > 15 GeV for events with MIPs. 

The e+e" — > qq(7) process is a potentially large source of background. It is reduced by 
exploiting the fact that it originates muons close to the jets and that the total missing mo- 
mentum, if any, points towards the beam direction. The product of ■^^j, the angle between the 
muon and the closest jet, and sin6'miss is required to be greater than 5.5 degrees for events with 
muons and greater than 20 degrees for events with MIPs. 

Background due to qqyU+/i~ final states from Z-boson pair production in events containing 
MIPs is rejected by requiring the relativistic velocity of the reconstructed W bosons to be 
greater than a y^-dependent value, ranging from 0.34 to 0.49. 

The purity of the selection is 98% at ^/s — 206.5 GeV. The residual background, not 
originating from W-boson pair production, is dominated by Z-boson pair-production events 
(52%) and e+e~ — >• qq(7) events (31%). The distributions of Mjj and of tjj^j x sin 6*111188 are 
shown in Figure 3. 

3.4 The qqriy Selection 

The event selection for the process e+e^ — > qqru is based on the identification of an isolated 
low-momentum electron, muon, or narrow jet in a hadronic environment with large missing 
energy. 

Events are selected requiring Pt > 10 GeV, 30 GeV < Mjj < 110 GeV and the mass recoiling 
against the two-jet system to be greater than 35 GeV. 

Events are classified according to the presence of isolated electrons or muons with an energy 
of more than 5 GeV. MIPs are not considered as tau candidates. 

For leptonically-decaying tau candidates, cuts on Mei, and M^,^, complementary to those de- 
scribed in Sections 3.2 and 3.3, are applied. These cuts are chosen so as to minimise correlations 
in the measured W-boson branching fractions. 

If no electrons or muons are found, a search for a tau-jet is performed using a neural network 
which exploits the distinctive characteristics of a hadronic tau decay: low multiplicity, small jet 
opening angle, low jet mass and high electromagnetic fraction of the jet energy. The jet with 
the highest neural-network output is retained as the tau-lepton candidate. For these events, 
additional requirements are applied in order to reduce the dominant background from e"'"e~ — >■ 
qq(7) events. If Pt < 20 GeV, the neural-network output of the tau-jet candidate is required 
to be near to that expected for a tau-jet. At most three charged tracks are allowed to form the 
tau-jet candidate. The polar angle of the missing momentum must satisfy | cos^missl < 0.91. 
The solid angle subtended by the directions of the tau-jet candidate and the other two jets 
must be less than 6 sr. 

Among the events selected at ^/s = 206.5 GeV, 62% come from e"''e~ — > qqru W-boson 
pair-production processes and 21% from other final states of the W-boson pair production. 
The background is dominated by e+e~ — > qq(7) events (54%) and e+e" — > qqeu final states not 
originating from W-boson pair production (46%) . The distributions of M^,^ and Mjj are shown 
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in Figure 3. 



3.5 The qqqq Selection 

The event selection for the process e''"e~ — > qqqq requires hadronic events with httle missing 
energy, high multiphcity and a four-jet topology. 

The Durham jet-resolution parameter 7/34. for which the event topology changes from three 
to four jets, is required to be greater than 0.0015. The events are clustered into four jets and 
a kinematic fit, assuming four-momentum conservation, is used to improve energy and angle 
resolutions. 

A neural network is trained to discriminate against the dominant e"'"e~ — > qq(7) background. 
Ten variables are used in the neural network: the spherocity [34], the lowest jet-multiplicity, 
1/34, the energies of the most and of the least energetic jets, the difference between the energies 
of the second and the third most energetic jets, the broadenings [35] of the most and of the 
least energetic jets, the probability of the kinematic fit and the sum of the cosines of the six 
angles between the four jets. The dominant background is due to e"''e~ — qq(7) events with 
four reconstructed jets, mainly coming from e+e" qqgg events. We find that the four-jet 
rate in e"'"e~ qq(7) events is not well described by our MC simulations, and a comparison 
with data is used to determine this background. Data and Monte Carlo distributions of the 1/34 
variable in hadronic Z decays collected at ^/s — 91 GeV are compared and their ratio is used to 
reweight the e"'"e~ — > qq(7) Monte Carlo events at higher energies throughout the rest of the 
analysis. The resulting accepted number of e"'"e~ qq(7) events, for a neural network output 
greater than 0.6, is increased by 12.7%. 

Requiring the neural-network output to be greater than 0.6 yields a sample purity of 80% 
at ^/s — 206.5 GeV with a background dominated by the e+e~ qq(7) (59%) and Z-boson 
pair-production (41%) processes. The distributions of some of the neural-network inputs and of 
its output, peaking at one for the signal and at zero for the background, are shown in Figure 4. 



4 Fit Method 

The CC03-level cross sections, aj, of the signal processes j are determined simultaneously in a 
single maximum- likelihood fit, taking cross-feed between different final states into account. 

For the purely leptonic final states, the fit procedure determines six different cross sections 
corresponding to all possible lepton-fiavour combinations. Since the statistics for the £i/£i/ final 
state is low, the sum of these six cross sections is quoted in the following as the cross section 
for the process e"'"e~ — > iuiu. 

The total likelihood is given by the product of Poissonian probabilities, P{Ni, fj^i), to observe 
Ni events in the i-th final state, as listed in Table 3. The expected number of events for selection 
i, /li, is calculated as: 

/ 10 "^^'b b ^ 

\j=i k=i J 

where Sij is the CC03-level efficiency of selection i to accept events from process j, a^^ is the 
cross section of the k-th background process, selected with efficiency e^f. The ^ background 
processes for selection i also include four-fermion final states not originating from W-boson pair 
production. 
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For the e+e" qqqq process, the Poissonian probabihty is replaced by the hkehhood as a 
function of the signal cross section derived from a fit to the neural-network output distribution. 
In this fit the e+e^ qq(7) background contribution is fixed to the value derived directly from 
data by performing a fit with both the signal and background normalisations left free. The 
results for the e+e~ qq(7) background cross sections are shown in Table 4. These values are 
in good agreement with the Monte Carlo predictions. As a cross-check, the e+e~ — > qqqq cross 
section is also determined by repeating the full fit after applying a cut on the output of the 
neural network at 0.6, which minimises the expected statistical uncertainty. All values agree 
well with those derived from the neural-network fit. 

5 Systematic Uncertainties 

In addition to the uncertainty on the luminosity measurement [11] and that due to limited Monte 
Carlo statistics, which affect all final states in common, the remaining sources of systematic 
effects in the measurement of W-boson pair-production cross sections are divided into two 
classes: uncertainties in the detector response and theoretical uncertainties. The latter come 
mainly from the knowledge and modelling of the hadronisation processes. A summary of the 
systematic uncertainties from all considered sources is given in Table 5 for ^/s — 206.5 GeV. 
Values at different ^/s are only marginally different. Details about the assessment of the 
systematic uncertainties are discussed below. 

A possible source of systematic uncertainty arises from the accuracy of the Monte Carlo 
modelling of the detector response. For the semi-leptonic and fuUy-leptonic final states, this 
uncertainty is evaluated by varying the positions of the selection cuts for each channel. The 
variation of the cut positions is chosen so as to span several times the resolution of the studied 
variable. Each variable is considered in turn and the corresponding change in the measured 
cross sections are evaluated. For variables which are correlated, for instance visible energies and 
transverse momenta, the largest variation is retained. For the selected variables, the expected 
statistical uncertainty on the newly-selected data sample is subtracted from the observed vari- 
ation and the sum in quadrature of all results is retained as systematic uncertainty. Most of 
the systematic effects are related to the resolution of the missing momentum. In addition, the 
electron/photon discrimination represents also a sizable source of systematic uncertainty for 
the qqeu final state. 

For the qqqq selection, the systematic uncertainty on the neural-network output is esti- 
mated by re-evaluating the input variables of the neural network after smearing and scaling the 
measurements of energy depositions and tracks in the simulation according to the uncertainties 
on their resolutions. 

The relative systematic uncertainty on the measured cross section, assigned to detector 
response and modelling, varies from 1.0% to 2.0% depending on the final state. 

As a cross-check, changes in efficiency due to variations of the detector calibration within its 
uncertainty, are also studied. The calibration is studied using samples of di-lepton and di-jet 
events, collected during the calibration runs at i/i = 91 GeV and at higher energies. The 
results of this study show a much smaller effect than the cut-variation technique. The trigger 
inefficiency, as well as its uncertainty, is found to be negligible in all channels. 

Fragmentation and hadronisation uncertainties may affect both the signal efficiency and the 
e"'"e~ — > qq(7) background estimation. The modelling of the signal hadronisation is studied 
comparing the selection efficiencies obtained with different hadronisation models: PYTHIA, 
HERWIG [36] and ARIADNE [37]. The average difference with respect to PYTHIA gives a 
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systematic uncertainty on the measured cross section of 0.5% to 1.2%, dependent on the final 
state. 

The effect of the hadronisation uncertainty in e+e~ '¥l{l) background events is also 
studied by comparing PYTHIA, HERWIG and ARIADNE. It is found to be negligible for qq£i/ 
final states. In the qqqq final state, the hadronisation uncertainty affects mainly the four-jet 
rate as described in Section 3.5. Half of the effect due to the reweighting is assigned as 
systematic uncertainty. It corresponds to 0.9 % of the measured c+c" qqqq cross section. 

Other sources of theoretical uncertainties in the qqqq channel arise from correlations among 
final-state hadrons such as Bose-Einstein correlations and colour reconnection. The modelling 
of Bose-Einstein correlations between hadrons from W-boson decays may affect the selection 
efficiencies. In previous studies [29] we have measured the strength of Bose-Einstein correlations 
between hadrons originating from the same W boson in semi-lcptonic W decays. Its value is 
significantly different from zero and in good agreement with that for light-quark Z decays and 
also with that of the LUBOEI BE32 model [28] used in our Monte Carlo simulations. The 
systematic uncertainty derived from the uncertainty of this strength is found to be negligible. 
Bose-Einstein correlations between particles originating from different W bosons are strongly 
disfavoured in e+e" qqqq events [29]. Their measured strength is restricted to at most a 
quarter of the strength expected in the BE32 model with full correlations. Allowing correlations 
with such a strength yields negligible changes in the measured cross sections. 

Extreme models of colour reconnection between the hadronic systems in qqqq events are 
disfavoured by data [38, 39] . The infiuence of colour reconnection is estimated using the models 
implemented in HERWIG [40], ARIADNE [41] (model 1 and model 2) and PYTHIA (model SK 
I with reconnection parameter k = 0.6 [42]). The ARIADNE-2 model is compared to a modified 
version of the ARIADNE- 1 model, so that in both models the shower cascade is performed in 
two phases with an identical energy cut-off parameter. The average difference of 0.19% with 
respect to PYTHIA is assigned as systematic uncertainty on the measured e+e~ — > qqqq cross 
section. 

The theoretical uncertainties on the cross sections of the background processes, namely 
hadron production in two-photon collisions (50%), neutral-current four-fermion processes (5%) 
and fermion-pair production (1%) lead to systematic uncertainties of 0.1% to 0.4%. In the 
determination of the e"'"e~ — > qqqq cross sections, the e"'"e~ — > qq(7) background levels are 
directly measured from data and the corresponding uncertainties, as reported in Table 4, are 
propagated to the final results. 

The dependence of the selection efficiencies on the mass and width of the W boson, mw 
and Fw, is studied using Monte Carlo samples simulated with different mw and Fw values. 
The propagation of the world-average uncertainties on these two parameters, 40 MeV on mw 
and 60 MeV on Fw [43], is taken as systematic uncertainty. It corresponds to a less than 0.3% 
effect. 

The systematic uncertainty on initial-state radiation (ISR), due to its approximate leading- 
log 0{a^) treatment in KandY, is investigated by re-evaluating the signal efficiencies for Monte 
Carlo events reweighted by |A^[0(q;^)]P/|A^[0(q;^)]|^. The effect is found to be neghgible. 
As a cross-check, the Monte Carlo events are also reweighted by 10% in the presence of ISR 
photons with energies or transverse momenta exceeding 100 MeV. In both cases, the effect is 
negligible. 

Final-state radiation (FSR) is implemented in KandY using the PHOTOS package [44] 
based on the leading-log approximation. The PHOTOS package is inaccurate in the hard non- 
coUinear region. The related systematic uncertainty is estimated by determining the selection 
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efficiencies using Monte Carlo events wliose weiglits are reduced by 50% in tlie presence of FSR 
pliotons witli energy greater tlian 30 GeV. An effect between 0.1% and 0.2% is observed and 
retained as systematic uncertainty. 

Uncertainties due to tlie implementation of virtual 0{a) corrections in the KandY program 
are tested comparing signal efficiencies to those obtained with the Racoon WW program. No 
sizable effect is observed. 

Correlations among all sources of systematic uncertainties are taken into account in the 
following results. 

6 Results 

6.1 Single- Channel Cross Sections 

Fits are performed to derive ten cross sections, one for each final state. No assumption is made 
concerning the W-boson branching fractions. The results, summing up all fully leptonic final 
states and including statistical and systematic uncertainties, are listed in Table 3. The Standard 
Model agrees well with these results. Since the efficiency matrix of Table 2 contains non-zero off- 
diagonal elements, the measured cross sections are correlated. The largest correlations, —10.3% 
and —17.6%, are between the e+e" qqru and e+e~ — >• qqeu and between the e+e" — > qqri/ 
and e+e~ — > qq/iu cross sections, respectively. All other correlations are less than 1%. 

6.2 Total Cross Section and Branching Fractions 

For the determination of the CC03 cross section of W-boson pair production, o"ww, the signal 
cross sections aj are replaced by the product r^aww- The ratios rj are given in terms of the W- 
boson decay branching fractions, Br(W — > qq) and Br(W iu), as follows: rqqqq = [Br(W 
qq)]^, r^qii, — 2Br(W — > qq)Br(W — > iu), and vmi, — [Br(W iv)Y for same-flavour leptons 
or 2Br(W £z/)Br(W ^ f z/) otherwise. 

Results for the cross sections of the reactions c+e~ — > iviv^ e+e" qqiv and e+e" qqqq, 
assuming charged- lepton universality, are obtained as shown in Table 6. The total cross sections, 
(Tww; are then derived assuming the Standard Model W-boson decay branching fractions [3] 
and are also reported in Table 6 together with the Standard Model expectations. Our previous 
measurements at y^s of 161 GeV [5], 172 GeV [6], 183 GeV [7] and these results are compared 
in Figure 5 to the Standard Model expectation as calculated with the Monte Carlo programs 
YFSWW3 and Racoon WW. The two predictions agree with our data and are consistent within 
a theoretical uncertainty of 0.5% [45] for ^/s > 170 GeV. 

The ratios of the measured cross sections to the Standard Model predictions of the YFSWW3 
program are also shown in Figure 5. Their combined value, R, is: 

i?= 0.992 ±0.011 ±0.009 ±0.005, (2) 

where the first uncertainty is statistical, the second systematic and the third theoretical. 

For the determination of W-boson decay branching fractions, the data collected at lower 
centre-of-mass energies are also included. The sum of the hadronic and the three leptonic 
branching fractions is constrained to unity. The branching fractions are first determined without 
the assumption of charged-lepton universahty, with the results hsted in Table 7. The hypothesis 
of charged-lepton universality is tested and the probability of getting a greater than that 
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observed is 0.8% differing by 2.6 standard deviations from tliis fiypotfiesis. Assuming cfiarged- 
lepton universality, the fiadronic W-boson decay branching fraction is: 

Br(W^qq) = 67.50 ± 0.42 ± 0.30 % , (3) 

where the first uncertainty is statistical and the second systematic. The W-boson decay branch- 
ing fractions depend on the six elements Vij of the Cabibbo-Kobayashi-Maskawa matrix V [46] 
not involving the top quark [3] : 

l/Br(W^£i/) = 3 + 3[1 + a,(mw)/7r] J2 l^^il'' 

i=u,c; j=d,s,b 

where ccg is the strong couphng constant. Using = 0.119 ± 0.002 [43], our measurements 
correspond to: 

J2 \Vijf = 2.002 ±0.038 ±0.027, (4) 

i=u,c; j=d,s,b 

where the first uncertainty is statistical and the second systematic. 

Using the current world-average values and uncertainties of the other matrix elements, not 
imposing the unitarity of the V matrix, |V^s| is derived as: 

iKsl = 0.977 ±0.020 ±0.014, (5) 

where the first uncertainty is statistical and the second systematic. The systematic uncertainty 
includes the uncertainties on ctg and on the other matrix elements [43]. 

6.3 Differential Cross Section 

The combined differential cross section for the e+e^ qqez/ and e+e~ — > qq/xi/ channels, as a 
function of cos ^w- > where 9^^- is the W~ production angle with respect to the direction of the 
incoming electrons, is measured for different ^/s from 183 GeV to 209 GeV. These two channels 
arc used because the lepton charge tags the W-boson charge with high purity. 
Four energy bins are considered: 

180.0 - 184.0 GeV, 184.0 - 194.0 GeV 
194.0 - 204.0 GeV, 204.0 - 209.0 GeV . 

These are chosen so as to minimise the difference between the average slope of the differential 
cross section in each bin and the slope corresponding to the luminosity-weighted average centre- 
of-mass energies: {^/s) = 182.7, 189.0, 198.3 and 205.9 GeV, respectively. In each energy range, 
ten cos^w- bins are studied. The variable cos^w- is reconstructed from the measurements 
of the jet and lepton angles and energies [47]. Monte Carlo events are then used to extract 
the differential cross section. Ambiguities might arise in the presence of additional photons 
in the generated events, and the cos^w- angle is then defined following the 7-recombination 
scheme [45]: 

• all photons inside a cone of 5 degrees half-opening angle with respect to the beam direction 
are treated as invisible; 

• the combined mass of each photon with electrons, muons and quarks is calculated. If the 

smallest combined mass is less than 5 GeV or the energy of the photon is less than 1 GeV, 
the momentum of the photon is added to that of the fermion and the photon is discarded. 
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The measured cross sections are corrected to CC03-level with the additional restriction of 
20° < 9^± < 160°, where 6'^± is the angle between the charged lepton and the beam direction. 

The observed cos^^w- distributions are corrected to generator level, after background sub- 
traction, by using bin-by-bin correction factors and the cross sections in each cos6'w- bin are 
determined as hsted in Table 8 and plotted in Figure 6. 

As a cross-check, a full matrix unfolding from reconstruction to generator level is also used. 
Since the migration matrix is almost diagonal, with bin-to-bin migration effects at the level of 
20% at most, the results are in perfect agreement with the simple bin-by-bin correction method. 

The potential bias of implicitly assuming the Standard Model cos 9^^- distribution in the 
correction factors, is studied using simulated samples with modified cos^w- behaviour and 
found to be negligible. Another bias could arise directly from the W-boson pair-production 
Monte Carlo generator used to estimate the correction factors. No difference between KORALW 
and YFSWW3 programs is observed, hence no additional systematic uncertainty is assigned. 

Charge-confusion effects, which affect the reconstruction of the W-boson direction, are taken 
into account. The residual uncertainty, obtained by comparing data and Monte Carlo expec- 
tations on Z-peak samples [47], is retained as a systematic uncertainty in addition to those 
affecting the total W-boson pair-production cross section. 

The systematic uncertainty is taken to be fully correlated between cos 9^^- bins and energy 
points. 

7 Conclusions 

In a data sample corresponding to an integrated luminosity of 629.4 pb~^, collected at centre- 
of-mass energies ranging from 189 GeV to 209 GeV, W-boson pair-production cross sections are 
measured by selecting four-fermion events and found to be in agreement with Standard Model 
expectations. 

The branching fractions for leptonic W decays are measured for each lepton generation. 
Assuming charged-lepton universality, the branching fraction for hadronic W decays is measured 
to be: 67.50±0.42 (stat.)ib0.30 (syst.) %. Combining all ^/s points, the ratio R of the measured 
total W-boson pair-production cross section with respect to the theoretical prediction is found 
to be i? = 0.992 ± 0.011(stat.) ± 0.009(syst.) ± 0.005(thco.). 

Differential cross sections as a function of the W^ production angle are also measured and 
found to be in good agreement with Standard Model predictions. 
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(v^) [GcV] 


188.6 191.6 195.5 199.6 201.8 204.8 206.5 208.0 




176.8 29.8 84.1 83.3 37.1 79.0 130.5 8.6 



Table 1: Average centre-of-mass energies and integrated luminosities. 



Selection 


Efficiencies [%] for e^ 
eueu eufiu euTU iiuiiu fiuTU tutu 


^e — * 
qqeu qq/iu qqTU 


qqqq 


evev 

TUTU 


54.7 0.8 11.4 0.1 1.5 
47.6 8.4 1.4 10.1 2.2 
6.0 1.7 27.8 0.4 7.5 
41.0 6.9 0.9 
2.6 0.3 3.0 23.1 4.8 
0.2 0.1 2.1 1.3 16.7 






qqcu 
qqfiu 
qqTU 




73.3 0.2 1.6 
0.1 74.2 4.2 
6.2 10.1 49.8 


0.1 


qqqq 




0.1 0.4 


84.0 



Table 2: Selection efficiencies for the signal processes e"'"e~ — > iuiu, e'^e~ qqiu, and e+e" — > 
qqqq, at ^/s = 206.5 GeV. For the e"'"e~ — > qqqq selection, the numbers are quoted for a 
neural- network output greater than 0.6. Selection efficiencies at other centre-of-mass energies 
are only marginally different. 
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[pb] 


[pb] 




(Vs) = 188.6 GeV 


(Vs) = 201.8 GeV 


iviv 


235 


57.2 


1.87 ±0.17 ±0.06 


1.72 


40 


12.3 


1.47 ±0.35 ±0.07 


1.81 


qqei/ 


347 


22.9 


2.29 ±0.14 ±0.03 


2.38 


70 


5.3 


2.26 ±0.30 ±0.03 


2.49 


qq//i/ 
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14.9 


2.25 ±0.14 ±0.04 


2.38 


79 


3.4 


2.62 ±0.33 ±0.05 


2.49 


qqri/ 


413 


69.7 


2.82 ± 0.22 ± 0.07 


2.38 


77 


13.9 


2.45 ±0.47 ±0.06 


2.49 


qqqq 


1477 


328.7 


7.17 ±0.24 ±0.12 


7.42 


301 


64.6 


7.10 ±0.52 ±0.12 


7.79 




(Vs) = 191.6 GeV 


(Vi) = 204.8 GeV 


iviv 


35 


10.4 


1.67 ±0.41 ±0.07 


1.76 


85 


25.9 


1.58 ±0.26 ±0.05 


1.82 


qqei/ 


73 


4.1 


2.95 ±0.37 ±0.04 


2.42 


176 


11.0 


2.78 ±0.23 ±0.04 


2.50 


qq//i/ 


63 


2.4 


2.61 ±0.36 ±0.04 


2.42 


142 


6.5 


2.30 ±0.22 ±0.04 


2.50 


qqri/ 


57 


11.9 


1.87 ±0.48 ±0.05 


2.42 


164 


26.4 


2.63 ±0.33 ±0.07 


2.50 


qqqq 


236 


57.5 


6.79 ±0.56 ±0.15 


7.56 


656 


137.2 


7.66 ±0.37 ±0.13 


7.81 




(v^) = 195.5 GoV 


(v^) = 206.5 GeV 


iviv 


105 


30.2 


1.76 ±0.25 ±0.06 


1.79 


128 


42.6 


1.42 ±0.19 ±0.06 


1.82 


qqei/ 


168 


10.9 


2.36 ±0.20 ±0.03 


2.46 


269 


16.9 


2.56 ±0.17 ±0.03 


2.50 


qq//i/ 


157 


8.2 


2.14 ±0.20 ±0.03 


2.46 


240 


11.8 


2.28 ±0.17 ±0.04 


2.50 


qqri/ 


222 


33.8 


3.44 ± 0.34 ± 0.08 


2.46 


287 


45.1 


2.92 ±0.27 ±0.07 


2.50 


qqqq 


665 


153.5 


6.92 ±0.34 ±0.11 


7.68 


1108 


220.1 


8.07 ±0.29 ±0.13 


7.82 




(Vi) = 199.6 GeV 


(v/i) = 208.0 GeV 


iviv 


87 


26.0 


1.68 ±0.27 ±0.06 


1.80 


11 


2.4 


2.23 ±0.86 ±0.06 


1.82 


qqei/ 


152 


11.4 


2.21 ±0.20 ±0.03 


2.48 


14 


1.1 


2.02 ±0.61 ±0.03 


2.50 


qq//i/ 


142 


7.3 


2.05 ±0.20 ±0.04 


2.48 


23 


0.7 


3.59 ±0.81 ±0.05 


2.50 


qqri/ 


181 


32.2 


2.75 ± 0.32 ± 0.07 


2.48 


17 


2.9 


2.43 ± 1.03 ±0.06 


2.50 


qqqq 


726 


151.1 


7.91 ±0.36 ±0.13 


7.76 


65 


14.1 


7.28 ± 1.16 ±0.11 


7.82 



Table 3: Number of selected data events, A^data, number of expected background events, A^bg, not 
originating from W-boson pair production, and CC03 cross sections for the reactions e'''e~ — > 
Ivlv^ e"'"e~ qqez/, e+e" — * qqAiz/, e+e" — > qqrz/, and e+e" qqqq. For e+e~ — > qqqq, 
-^data and iVbg correspond to a cut on the output of the neural network at 0.6, while the 
e+e" qqqq cross section is obtained from a fit to the neural-network output distribution, as 
described in Section 4. All cross sections are derived without any assumption on the W-boson 
decay branching fractions. The first uncertainty is statistical and the second systematic. Also 
shown are the Standard Model CC03 cross sections, usm, as calculated with YFSWW3 [16] 
with an uncertainty of 0.5%. 
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(v/i) [GeV] 


a--^ [pb] 


<^ [pb] 


188.6 


107.5 ±3.4 


101.00 


191.6 


92.8 ±7.6 


97.74 


195.5 


86.7 ±4.5 


92.47 


199.6 


86.8 ±4.7 


88.09 


201.8 


89.6 ±7.0 


85.89 


204.8 


84.1 ±4.7 


82.19 


206.5 


78.1 ±3.6 


80.90 



Table 4: Measured, cr™^^, and expected, cr^, cross sections of the e"''e" — > qq(7) process. The 
measurements are determined by a fit of the neural- network output distribution of the qqqq 
selection to both signal and e+e~ —>■ qq(7) background. 



Systematic uncertainties on a [% 



Final state 



Source 




qqev 




qqrv 


qqqq 


Luminosity 


0.22 


MC statistics (signal) 


0.80 


0.25 


0.25 


0.44 


0.11 


MC statistics (background) 


1.57 


0.23 


0.28 


0.75 


0.22 


Detector modelling 


2.00 


1.00 


1.20 


2.00 


1.00 


Hadronisation (signal) 




0.77 


0.58 


1.17 


0.45 


Hadronisation (background) 










0.90 


Bose- Einstein effects 




< 0.01 


0.03 


Colour reconnection 










0.19 


Background cross sections 


0.59 


0.21 


0.22 


0.40 


0.40 


W mass (±0.04 GeV) 


0.27 


0.03 


0.03 


0.10 


0.06 


W width (±0.06 GeV) 


0.12 


0.03 


0.12 


0.08 


0.02 


ISR simulation 


< 0.01 


FSR simulation 


0.21 


0.21 


0.17 


0.08 


< 0.01 


Total 


2.76 


1.36 


1.43 


2.52 


1.46 



Table 5: Relative systematic uncertainties on the cross-section measurements evaluated for 
a/s — 206.5 GeV. Uncertainties at other center-of-mass energies are only marginally different. 
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r 1 1 

a [pbj 


(Vs) = 188.6 GcV 


{^ys) = 191.6 GcV 


{y/s) = 195.5 GcV 


/ 1 \ 1 r\r\ i r 

(v's) = 199.6 GcV 


O'qqqq 


1.00 ± U.16 ± (J. (J / 

7.19 ±0.24 ±0.08 
7.17 ±0.24 ±0.12 


1.66 ± U.o9 ± U.U 1 

7.69 ±0.61 ±0.09 
6.78 ±0.56 ±0.12 


1. /o ± U.24 ± U.U 1 

7.58 ±0.36 ±0.08 
6.92 ±0.34 ±0.11 


1. ± U.ZD ± U.U6 

6.81 ±0.35 ±0.08 
7.91 ±0.36 ±0.13 


0"WW 


1 ^ 1 T _i_ n o v _i_ n 1 v 
lO.l / ± U.o ( ±[}.i( 


Id. 11 ± U.oy ± U.l / 


10.22 ± U.o4 ± U.lO 


lb. 49 ± U.OO ± U.l/ 


CSM 


lb. 2 1 


1 ^ rr V 

lo.o 


1 ^ O /I 


1 7 no 
1 i .\j2 


Til 

a [pb] 


[y/s) = 201.8 GcV 


(y/s) = 204.8 GcV 


(Vs) = 206.5 GcV 


(Vs) = 208.0 GcV 


^qqqq 


1.51 ±0.34 ±0.07 
7.34 ±0.54 ±0.08 
7.09 ±0.52 ±0.12 


1.58 ±0.24 ±0.05 
7.68 ±0.39 ±0.13 
7.66 ±0.37 ±0.13 


1.44 ±0.18 ±0.06 
7.60 ±0.30 ±0.08 
8.07 ±0.29 ±0.13 


2.23 ±0.86 ±0.06 
8.18 ± 1.21 ±0.09 
7.29 ± 1.16 ±0.11 


0"WW 


16.01 ±0.81 ±0.17 


17.00 ±0.58 ±0.17 


17.31 ±0.45 ±0.18 


17.52 ± 1.81 ±0.17 


CSM 


17.08 


17.12 


17.14 


17.15 



Tabic 6: Measured CG03 cross sections of the processes e+e" iviv^ e+e" — > qq£z/ (summed 
over lepton flavours) and e+e^ — > qqqq, assuming charged-lepton universality. The measured 
W-boson pair-production cross sections, crww, are derived assuming Standard Model branching 
fractions for the W boson decay modes. The Standard Model total W-boson pair-production 
cross sections, (Tsm, are calculated using the YFSWW3 program, which has a theoretical un- 
certainty of 0.5%. 



Branching 
fraction 


Lepton 
non- universality 


Lepton 
universality 


Standard 
Model 


Br(W ^ Qv) [%] 
Br(W ^ [iv) [%] 
Br(W ^ rz/) [%] 
Br(W ^ tv) [%] 
Br(W ^ qq) [%] 


10.78 ±0.29 ±0.13 
10.03 ±0.29 ±0.12 
11.89 ±0.40 ±0.20 

67.30 ± 0.42 ± 0.30 


10.83 ±0.14 ±0.10 
67.50 ± 0.42 ± 0.30 


10.83 
67.51 



Table 7: W-boson decay branching fractions derived without and with the assumption of 
charged-lepton universality. The correlation coefficients between the leptonic branching frac- 
tions are -0.016, -0.279, -0.295 for [Br(W ^ ez/), Br(W ^ \xv)\ [Br(W ^ ez/), Br(W ^ rz/)] 
and [Br(W — > yU.z/),Br(W — > rz/)], respectively. The W-boson decay branching fractions ex- 
pected in the Standard Model are also listed. 
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dcr/dcos^w- [pb] 


COS ^w- 


range 


/ /TTX 1 QO 7 

VV-^/ — i-oZ.( 


Qi\/r 


[y/S) — ioy.U 




-l.U, 


— U.o 


U.o4 ± U.zo ± U.Ui 


n 7/1 
U./4 


n Rn _L n 1 o _L n ni 


U.o4 


u.o, 


— n fi 
u.u 


fi Qi _|_ fi on _|_ fi ni 




nQQ_|_nic:_|_nni 




-0.6, 


-0.4 


0.22 ±0.26 ±0.00 


1.02 


1.08 ±0.17 ±0.02 


0.94 


-0.4, 


-0.2 


0.96 ±0.33 ±0.01 


1.20 


1.18 ±0.19 ±0.02 


1.14 


-0.2, 


0.0 


1.71 ±0.43 ±0.03 


1.44 


1.34 ±0.20 ±0.02 


1.38 


0.0, 


0.2 


2.27 ±0.50 ±0.03 


1.78 


1.51 ±0.22 ±0.02 


1.72 


0.2, 


0.4 


3.37 ±0.62 ±0.05 


2.16 


1.88 ±0.24 ±0.03 


2.22 


0.4, 


0.6 


3.52 ±0.66 ±0.05 


2.86 


2.95 ±0.31 ±0.04 


2.95 


0.6, 


0.8 


4.24 ±0.74 ±0.06 


3.84 


4.19 ±0.37 ±0.06 


4.15 


0.8, 


1.0 


5.00 ±0.83 ±0.07 


5.47 


6.11 ±0.47 ±0.09 


6.24 





dcr/d cos 


^w- [pb] 




cos ^^w- 


range 


(^) = 198.3 


SM 


(v^) = 205.9 


SM 


-1.0, 


-0.8 


0.68 ±0.11 ±0.01 


0.57 


0.60 ±0.10 ±0.01 


0.52 


-0.8, 


-0.6 


0.76 ±0.13 ±0.01 


0.71 


0.44 ±0.11 ±0.01 


0.64 


-0.6, 


-0.4 


0.78 ±0.15 ±0.01 


0.85 


0.77 ±0.14 ±0.01 


0.78 


-0.4, 


-0.2 


0.80 ±0.16 ±0.01 


1.05 


0.99 ±0.16 ±0.01 


0.98 


-0.2, 


0.0 


1.31 ±0.20 ±0.02 


1.29 


1.35 ±0.20 ±0.02 


1.21 


0.0, 


0.2 


1.64 ±0.23 ±0.02 


1.65 


1.72 ±0.23 ±0.03 


1.55 


0.2, 


0.4 


2.21 ±0.27 ±0.03 


2.16 


1.75 ±0.23 ±0.03 


2.06 


0.4, 


0.6 


2.41 ±0.29 ±0.04 


2.97 


2.84 ±0.30 ±0.04 


2.92 


0.6, 


0.8 


3.69 ±0.36 ±0.05 


4.38 


4.80 ±0.41 ±0.07 


4.45 


0.8, 


1.0 


6.26 ±0.49 ±0.09 


7.20 


7.49 ±0.53 ±0.11 


7.80 



Table 8: Sum of the differential cross sections, as function of cos^w-) fo'^ the e+e^ — > qqeu 
and e'*"e~ — > qq/iiy processes. The measurements are derived in a restricted phase space of 
the CC03 subset of diagrams. The first uncertainty is statistical and the second is systematic. 
The systematic uncertainty is fully correlated between cos^^w- bins and between ^/s bins. 
The columns labeled SM show the expected values from the Standard Model, which have a 
theoretical uncertainty of about 2%. 
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e 




Figure 1: The lowest-order Feynman diagrams (CC03) contributing to W-boson pair produc- 
tion: t-channel v exchange and s-channel 7 and Z-boson exchange. 
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Figure 2: Distributions of variables used for the selection of iviv and qqei^ events, comparing 
the signal and background Monte Carlo to the data collected at -^i = 189 — 209 GeV. The 
positions of the selection cuts are indicated by vertical arrows. All selection cuts except the 
one on the plotted variable are applied, (a) The acoplanarity between the two leptons in 
the lepton-lepton class of the Ivlv selection. The excess in the first bin is due to cosmic-ray 
background, (b) The momentum transverse to the beam direction of the selected iviv events 
in the lepton-lepton class, (c) The energy of the electron in qqez/ events, (d) The absolute 
value of the cosine of the polar angle of the missing momentum in qqez/ events. 
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Figure 3: Distributions of variables used for the selection of qq/ii/ and qqrz/ events, comparing 
the signal and background Monte Carlo to the data collected at = 189 — 209 GeV. (a) The 
invariant mass of the jet-jet system in qq/iz/ events, (b) The quantity ip^j x sin O^^ss in qq/xi^ 
events, (c) The invariant mass of the jet-jet system in qqri/ events, (d) The invariant mass of 
the lepton-neutrino system for leptonically decaying tau candidates in qqrz/ events. 
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Figure 4: Distributions of some of the variables used for the neural network in the analysis of 
qqqq events together with the neural- network output, comparing the signal and background 
Monte Carlo to the data collected at = 189 - 209 GeV. (a) The spherocity. (b) The 
maximum jet energy scaled by the visible energy, (c) The broadening of the most energetic jet. 
(d) The neural- network output. 
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Figure 5: The cross section of the process e+e^ W^W~ as a function of ^/s. The pubhshed 
measurements of cxww at ^/s of 161 GeV, 172 GeV and 183 GeV, the updated measurement 
at ^/s = 189 GeV and the new measurements at ^/s = 192 — 209 GeV are shown as dots with 
error bars, combining statistical and systematic uncertainties in quadrature. The sohd curve 
shows the Standard Model expectation as calculated with YFSWW3 in the whole energy range 
and RacoonWW for -y/i > 170 GeV. Its uncertainty of 0.5% is invisible on this scale. The 
lower plot shows the ratios of the measured cross sections with respect to the Standard Model 
expectations as calculated with YFSWW3. The band represent their combined value with its 
total uncertainty: R = 0.992 ± 0.015. 
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Figure 6: Measured differential cross sections as a function of cos^w- for the e"'"e~ — > qqez/ and 
e"'"e~ —>■ qqfiu processes. Tfie cross sections of the two channels are summed. Experimental data 
are represented by dots with error bars which include statistical and systematic uncertainties 
added in quadrature. Monte Carlo expectations are shown as solid lines. 
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